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Abstract 

By making use of the double-time Green function technique, we study ther- 
modynamics of a deformed Bose gas, which describes well properties of density 
intensive photonic gas and radiation fields of the early universe. General form of 
statistical distribution function is obtained. We show explicitly the expression 
of the distribution function in some limitation cases. The free energy, equation 
of state, specific heat and other thermodynamic properties of the deformed Bose 
gas are presented. 



PACS numbers: 05.30.-d, 05.10.-a, 98.80.Cq, 05.70.-a, 51.30.+i 



1 Introduction 

There are two main motivations to study a deformed many-body system. The first one 
comes from fundamental physics. It is well-known that the noncommutative geometry 
plays a very important role in unravelling the properties of the Planck scale physics. 
It has for a long time been suspected that the noncommutative spacetime might be 

a realistic picture of how spacetime behaves near the Planck scale[l]. In fact, the 
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noncommutative geometry naturally enters the theory of open string in a background 
-B-field[2]. In particular, the noncommutative geometry makes the holography[3] {e.g. 
the AdS/CFT correspondence) of a higher dimensional quantum system of gravity and 
lower dimensional theory possible. A kind of special regularization with exponentially 
increasing spacetime cutoff has been introduced based on the noncommutative geome- 
try. A very small displacement of the noncommutative deformation parameter from its 
classical value reduces sharply the entropy of quantum system of gravity. An adequately 
adopted noncommutative deformation of geometry makes the holography of a higher 
dimensional quantum system of gravity and lower dimensional theory possible. It was 
also discovered that simple limits of M theory and superstring theory lead directly to 
the noncommutative gauge field theory [4, 5]. The noncommutative field theory has 
been intensively studied in the past two decads[6, 7]. The most realistic laboratory for 
testing the Planck scale physics is the early universe and the black hole physics. From 
the standard picture of cosmology, the universe started from a very hot, dense phase 
about 15 billion years ago in the big bang. The very early universe was opaque due to 
the constant interchange of energy between matter and radiation. Several problems of 
standard big bang cosmology, such as the horizon, homogeneity, entropy and flatness 
problems, can be explained under the assumption that the very early universe under- 
went a period of extremely rapid inflation[8], driven by the potential of some assumed 
inflation fields. In particular, the inflation scenario is also able to explain the observed 
perturbations in the cosmic microwave background radiation [9, 10], namely as origi- 
nating ultimately from quantum fluctuations of the inflation field. However, in most of 
the current models of inflation, the physical length of perturbations was much smaller 
than the Planck length at the beginning of inflationfll, 12]. It is then of interest to 
investigate whether the predictions of the spectrum of cosmological perturbations are 
sensitive to the unknown super- Planck- scale physics[13]-[16]. Recent studies[17]-[24] 
by physicists and astrophysicists show that without strange inflation fields but just 
the noncommutative geometry and the deformed radiation dominated universe subject 
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to the noncommutative spacetime quantization indeed give rise to inflation as the ra- 
diation temperature exceeds the Planck temperature. Based on the noncommutative 
geometry and the deformed radiation, one can even get an alternative explanation for 
the horizon, homogeneity, entropy and flatness problems in the standard cosmology. 

Another motivation comes from quantum optics as well as molecular and atomic 
spectroscopy. A fascinating development of quantum optics in last two decades has 
taken place within two specific areas, one concerning the study of nonclassical states 
of light, and the other studying the interaction of light with matter. The nonclassical 
states of light offer the possibility of the so-called reduction of quantum noise below the 
standard limit, determined by the Heisenberg's uncertainty relation. Study on interac- 
tion of light with matter shows that the old Bohr quantum jumps are real. Dynamics of 
deformed photonic field interacting with matter presented many interesting properties 
both in the two important areas [2 5] -[29]. Great success has also been made by making 
use of a deformed oscillator to study molecular spectrum[30]-[32]. 

In this paper, we study thermodynamics of a deformed Bose gas, and in particular, 
of a deformed radiation field. On the topic, many papers have been published in more 
than ten years [33]- [46]. However, we still have not a general form of statistical distri- 
bution function of the deformed Bose gas, which should reduce to the Bose-Einstein 
distribution function for an ideal Bose gas in the q — 1 limit. In this situation, it is dif- 
ficult to investigate systematically the influence of the deformation on thermodynamic 
characteristics of a Bose system, which is interesting in fundamental physics and other 
subjects as discussed above. For this purpose, by making use of the double-time Green 
function technique, we obtain an algebra equation, which should be satisfied by the 
distribution function of a deformed Bose gas. The equation for distribution function 
can not be solved exactly in general. The situation is similar to the case of magnetic 
lattices. From this point of view, a deformed Bose gas is really a strongly correlated 
system. In the case of q — > 1, we give a useful expression of the distribution function 
for a deformed Bose gas. Thermodynamic properties, such as internal energy, specific 
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heat, free energy and equation of state, are discussed in detail. 

The paper is organized as follows. In section 2, a deformed Bose gas is discussed. 
We get the distribution function for a deformed Bose gas by using the double-time 
Green function technique in section 3. In section 4, we calculate some thermodynamic 
quantities for a deformed photon gas and compare them with those for the ideal Bose 
gas. Concluding remarks are given in section 5. 

2 A deformed Bose gas system 

The Lagrange density of a scalar field 4>(x) with mass m reads 

h 2 d(fr d<p 1 2 2 2 

£{x) = — — — m c <p . 1 

v ; 2 dx^ dx» 2 r K J 

The Euler-Lagrange equation immediately leads to the Klein-Gordon equation 



The canonically conjugate field is 



d(p(x) 



which leads to the Hamilton density 

H{x) = \ (^tt 2 (x) + (W0(x)) 2 + m 2 c 2 2 (x)) . (4) 
The Fourier decomposition of the scalar field is of the form, 

(j)(x)= [ l d3k ° (a(k)e- ikx + aUk)e ikx ) , 

= I $k ih 2 u k ( a{k)e - ikx + a \ {k y kx \ 
J - /(27r) 3 2^ fe V ' 



where we have used the notation cu k = £yj (hk) 2 + m 2 c 2 . The canonical quantization 
of the scalar field reads 

[0(x,t),7r(y,t)]=^c5 3 (x-y) . (6) 



In terms of the Fourier modes a(k) and at(k), the commutation relations can be rewrit- 
ten into the form 

[a(k),at(k')] = 5 3 (k - k') . (7) 
It is not difficult to express the Hamiltonian in terms of these Fourier modes, 

H = -J Sk%u k (a{k)a\k) + a\k)a(k)) . (8) 

In the quantum field theory, the vacuum state is defined as follows 

a(Jfe)|0) = 0. (9) 

And the iV-particle Fock space can be written down directly as 

WkJMh), n(k m )) = v ^ v > y ■ ■ ■ v ) = |o) • (io) 

\jn[k x )\ \jn[k 2 )\ \jn[k m )\ 

To measure how many particles of a certain momentum exist in the field system, one 
can introduce the particle number operator N(k), 

N{k) = a\k)a(k) . (11) 

It is easy to check that 

N(k i )\n{k 1 ),n(k 2 ), ■ ■ ■ ,n(k m )) = n(k i )\n(k 1 ),n(k 2 ), ■ ■ ■ ,n(k m )) . (12) 

Thus, a field can be seen as a many-body system that consists of free harmonic oscilla- 
tors with different momentum. Properties of the field system can be described exactly 
by the Weyl-Heisenberg algebra, 

[a, a 1 "] = 1 , 

(13) 

[N, a] = -a , [N, a f ] = a) , 

with the Fock space \n) = ^=r|0) as its representation. All properties , in particular, 
the thermodynamic properties of the many-body system can be investigated system- 
atically, by making use of the Hamiltonian (8) and the Weyl-Heisenberg algebra (13). 

5 



The natural quantum of the quantum field theory on noncommutative geometry 
is the so-called deformed harmonic oscillator. On the topic of deformed oscillators, a 
lot of papers have been published[47]-[49]. Similar to the case of an ordinary quantum 
field system, the many-body properties of a noncommutative quantum field system can 
be described by a collection of the deformed harmonic oscillators, 

H q = \ E huJ k (a q (k)al(k) + a\(k)a q (k)) . (14) 

The deformed harmonic oscillator is related with the simple harmonic oscillator as 
follows [50], 

S a , a; = atS . (15) 



a q = 



q~ — q 

where we have used the notation [x] q = r- 

q — q 1 

The deformed Weyl-Heisenberg algebra reads 

a q a\ - qa\a q = q~ N , 

W, aj] = 4 ' t^' ^q] = ~a q ■ 
The representation of the deformed Weyl-Heisenberg algebra is obtained by construct- 
ing the Fock space, 

yW (17) 
a q \n) q = ^J\n\ q \n - l) q , a\\n) q = yj[n + l] q \n + 1), , 

where [n] q \ = [n] q [n — l] q ■ ■ • [2] q [l] q . It is easy to check that one still can measure how 

many particles of a certain momentum exist in the deformed field system by making 

use of the particle number operator N(k), 

N(k i )\n(k 1 ),n(k 2 ), ■ ■ ■ ,n(k m )) q = n{k i )\n(k 1 ),n(k 2 ), ■ ■ ■ ,n(k m )) q . (18) 

The Hamiltonian of a deformed harmonic oscillator [51] -[54] can be written as 

H q = ^([N] q + [N + l] q ). (19) 

Below we will calculate the distribution function of a deformed harmonic oscillator 
system by making use of the double-time Green function technique. 
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3 Distribution function 

The double-time Green function (retarded or advanced) for a Bose system [55, 56] is 
generally defined as 

G p (t - o =~0p+ l)0(t - + (p - l)0(f - 0) ( W), fl(01> - (20) 

where we have used the notations, 

i(*) = exp(if)iexp(=ip) , 

(6) = Z- lr Tr(e-P H d) , Z = Tr(e^ H ) , 

here p = 1 for the retarded one and p = — 1 for the advanced one. 
In Eq.(20), setting 

A(t) = a(t), s(0 = at(0, p = l, 

(the retarded Green function is assumed), we get 

G R (t - t') = -i0(t - t')(a(t)a\t') - a\t')a{t)) . (21) 
The equation of motion of the retarded Green function reads 
ih^-G R (t-t') = h (8{t - t'){a(t)a\t') - a\t')a(t)) + 9(t - t') (^^^(f) - a\t') 



da(t) 



dt 
(22) 



The Heisenberg equation for the operator a(t) 



gives, 

= -if + <T^ +1) ) a(t) . (23) 

Then, we have 



ih^G R (t - H) = fiS(t - t') + (q^ + q~^) G R (t - t') . (24) 



Making use of the Fourier transformation 

G R (t - = ^G R {E)e-^-^ h , (25) 

we get 

ft 

= ~ Mg^+g-^ • (26) 
The spectrum theorem of Green functions tells us 

(a (t )a(t)) - i J ^ dE _ i} . (27) 

Then, we can write formally the distribution function of a deformed Bose gas as 

n(e) = IN) = — 77T^ • (28) 

exp( ^^> +g -^>) ) _ 1 

Evidently, in the case of q — >1, n(e) recovers the familiar Bose-Einstein distribution 
function. In principle, all interesting properties of a deformed Bose gas can be ob- 
tained by solving this equation for distribution function. This is to be one of the most 
advantageous features of the Green function approach. 

Now we try to solve Eq.(28) to gain an insight into the influence of the deformation 
on the distribution function and to compare the deformed distribution function with 
the Bose-Einstein distribution function. 

Expanding the right hand side of Eq.(28) in the neighborhood of q = 1 to the second 
order, we get an exactly solvable equation for the deformed distribution function. The 
solution is of the form, 

+ 



5 2 f3eet 



^- £ \J(5 2 fiee^ + ( e /fe - l) 2 ) 2 - 5 2 (3ee^ {5 2 {3ee^ - 2{eP e - 1)) , 

(29) 

where we have used the notation 5 = q — 1 for short. 

We note that it is difficult to compare the distribution function of the deformed 
Bose gas, Eq.(29), with the Bose-Einstein distribution function directly. For the aim, 
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we rewrite Eq.(29) as follows, 

(e* _ 1) _ Wee* (d 2 (3ee^ - 2 (e* - l))' 



n e 



(30) 



(e* - l)^ + F/fee* (a^/fee* + 2 (e/3e _ 1)2 y 

One can further expand the distribution function into a series over deformation pa- 



rameter 8, 



n(e) 



1 



-8 



, , I3ee li£ I3ee' i£ pee 



/3e 



(31) 



V 2 ^" 1 ) 2 (e^-1) 3 2(e^-l)V ' 
From the above form of the distribution function of the deformed Bose gas, we see 
clearly that, at the case of q = 1, the Bose-Einstein distribution is recovered. To 
further investigate the effect of deformation on the distribution function, we present a 
plot of the distribution function in Fig.l. 




500 
400 
300 
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100 

0.01 0.02 0.03 0.04" 

FIG.l The distribution function of a deformed Bose gas 

( start from up q = 1.000, 1.001, 1.002, 1.003, 1.004, 1.005, 1.006, 1.010, here we have used the notation 0s = E.) 

To describe behaviors of the distribution function of the deformed Bose gas near the 
ground state more accurately, we give another plot of the deformed distribution function 
in Fig.2. 
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FIG. 2 The behaviors of distribution function 
of a deformed Bose gas near the ground state 

(start from up q = 1.000, 1.001, 1.002, 1.003, 1.004, 1.005, 1.006, 1.010, here we have used the notation fie = E.) 

4 Thermodynamic properties 

In this section, we concentrate on the study of thermodynamic properties of a deformed 
photonic field. As being well-known, the dispersion relation of a photon is 



e = cp . 



(32) 



We know that a photonic field can be described by a collection of simple harmonic 
oscillators and have relation 



H = hu[N + 



One can then read out that 



hu ( I 
p = — I n + 



(33) 



(34) 



c \ 1, 

This relation and the Hamiltonian (19) give the dispersion relation of a deformed Bose 
gas, 



e = -Tiuj 
2 



c 1 



+ 



hu P 2 



(35) 
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Expanding the right hand side of (35) in the neighborhood of q — 1 to the second 
order, we get 

S 2 5 2 , , 

£ = CP -24 CP+ W^ CP ■ (36) 
The density of energy levels is as follows, 

8nV SVWS + + 27Be 2 (-2 x 6^ A + (iSVBe + 2^3(4^ + 27 Be 2 )f") 



a(e) 



3h^ AB^AA* + 27W 2 ^ £ + 2y / 3{AA s + 27Be 2^ 1/3 

1/3 



9VBe + ^/3{AA 3 + 27Be 2 ) + &^A [18VBe + 2^3(AA 3 + 27Be 2 

X 2 

(9^^ + ^3(4^3 + 27Be 2 )) 

(37) 

where we have used the notations A = (1 -)c and .£> = — o-^-c 3 . 

V 24 ; 6ft V 

The mean value of the total energy of the system has the form, 



E 



/ en(e)a(e)de . (38) 
Jo 



To calculate E, we expand the right hand side of (37) to the third order in the neigh- 
borhood of q — 1 and rewrite the density of energy levels as the form 

Then we get 

where ((n) is the Riemann zeta function. 

We can further calculate the grand partition function In S through the thermody- 
namic relation E — — InS , 

op ' 

8n 5 V nVS 2 f8n 2 llvr 4 >/0 A /Im v, 

kS = " W (- + + «>C(3) J <^) 3 • (41) 

It is a straightforward calculation to obtain thermodynamic functions by making use 
of the grand partition function. The pressure of the deformed Bose gas is 
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The free energy of the deformed Bose gas is of the form, 



+ irV5 2 



six 



45 (he) 3 ' " ' " V 9 
The entropy of the deformed Bose gas is as 
U — F 



Uir 4 
~A5 



+ 



16C(3)1 



J (he) 3 



S 



3 ^ k4T3 v-^( 3 -^ + ^ + m s)) k ^v. 



(43) 



(44) 



T 45 h 3 c 3 V 9 45 7 

The first term in the above equation is just the entropy of an ordinary Bose gas. And 
the terms following it denote the effect of noncommutative deformation of geometry 
on the entropy. It is clear that the effect can not be ignored at least in the high 
temperature region, where the early universe lived. We give a plot of the entropy vs 
temperature for different deformation parameters in Fig. 3. 



5-10 



4 • 10 



3-10 



2 • 10 



1 • 10 




FIG. 3 The entropy per unit volume of a deformed Bose gas 

( start from up q = 1.0, 1.1, 1.2, 1.3, 1.4). 



Finally, the specific heat of the deformed Bose gas is 



32tt 5 V k A T 3 



'32n 2 Uir 4 
+ 



+ 192C(3) 



k 4 T 3 



(45) 



15 h 3 c 3 " ' " \ 3 ' 15 ' "~^"'J h 3 c 3 ' 
At the case of q — 1, we recover the result for an ordinary Bose gas. To show the effect 
of deformation on the specific heat, we give a plot of the specific heat per unit volume 
of a deformed Bose gas vs temperature for different deformation parameters in Fig. 4. 

From the above discusses, we can also obtain the Planck's formula for the deformed 
Bose gas, 
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2-10 4-10 6-10 8-10 1-10 
FIG. 4 The specific heat per unit volume of a deformed Bose gas 

( start from up q = 1.0, 1.1, 1.2, 1.3, 1.4). 




2 ' 4 ' 6 8 10 

FIG. 5 The Planck's formula of a deformed Bose gas 



( start from up q = 1.0, 1.1, 1.2, 1.3, 1.4, here we have used the notation / = f3Tu>j.) 

In Fig. 5, we give a plot of the Planck's formula for the deformed Bose gas. The first 
term in equation(46) describes an ideal gas consisting of photons. The linearity of 
the equations of electrodynamics on commutative spacetime guarantees that photons 
do not interact with one another. However, the other terms in this equation denote 
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interactions among photons or nonlinearity of equations of electrodynamics caused by 
noncommutative deformation of spacetime. At high frequencies {f3Twj ^> 1), we have 

hco 3 e-^ 175 2 hu; 3 e-P h " ^Ve"^ 

(47) 



u(uj, T) 



C 3 7T 2 



24c% 2 



2c% 2 



C 3 7T 2 2c 3 7T 2 

This is the Wien's formula of the deformed Bose gas. We give a plot of the Wien's 
formula for different deformation parameters in Fig. 6. In spite of a complex form of 
the Planck's formula and the Wien's formula, we note that the density of the spectral 
frequency distribution of the energy of the deformed Bose gas has a maximum at a 
frequency u m (f3hu m ~2.8), which is roughly independent of deformation parameter. 
Thus the displacement law is still remained for the deformed Bose gas. 




4 6 8 10 

FIG. 6 The Wien's formula of a deformed Bose gas 

( start from up q = 1.0, 1.1, 1.2, 1.3, 1.4, here we have used the notation / = (3hui). 



In the opposite limiting case of low frequencies ((3Tilo <C 1), we get the Rayleigh- Jeans 
formula for the deformed Bose gas, 



u(uj,T) = 



kTu 2 175 2 kTu 2 5 2 kTuj 2 ,hu 



2 24c% 2 

2(1 m\9 



2cH 2 y kT 



(77F + 1 ) 



5 2 {kT) 2 u huj_ ^ 5 2 {kTf (fiuj N 



(48) 



C 3 7T 2 /l ^kT 



2cH 2 h 2 \ kT 
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Fig. 7 shows a graph corresponding to the Rayleigh- Jeans formula with different defor- 
mation parameters. 

u 




f 



0.2 0.4 0.6 0.8 1 

FIG. 7 The Rayleigh-, Jeans formula of a deformed Bose gas 

( start from up q = 1.00, 1.06, 1.12, 1.18). 

We obtain the Stefan- Boltzmann's law for the deformed Bose gas of the form, 




0.2 0.4 0.6 0.8 1 



FIG. 8 The Stefan- Boltzmann law of a deformed Bose gas 

( start from up q = 1.0, 1.1, 1.2, 1.3, 1.4). 

Similar to the case of an ideal Bose gas, the total energy of the deformed Bose gas 
is proportional to the fourth power of the temperature. But the Stefan- Boltzmann's 
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constant is effectively reduced by the noncommutative deformation of geometry. The 
fact gives also a signal that the deformed Bose gas is a strongly correlated system. We 
give a plot of the Stefan- Boltzmann's law for different deformation parameters in Fig.8. 

In an adiabatic expansion (or compression) process, the entropy of the system is a 
constant and the valume and temperature are related by VT 3 =constant. From (42), 
we can get the equation of state, 

p\/ 4/3 = c q . (50) 

5 Concluding Remarks 

The noncommutative geometry plays a very important role in the study of the Planck 
scale physics. Quantum field theory and superstring theory (or M theory) on noncom- 
mutative spacetime have been studied for a long time. It is well-known that the most 
realistic laboratory for testing the Planck scale physics is the early universe. The very 
early universe was opaque due to constant interchange of energy between matter and 
radiation. It is clear that the radiation field at the early universe can not be described 
well by an ideal Bose gas system. There should be strong correlations among photons 
at this stage of evolution of the universe. In fact, attention has been paid recently to 
the investigation on the deformed radiation dominated universe. Preliminary results 
show that spacetime quantization gives rise to inflation as the radiation temperature 
exceeds the Planck temperature. Of course, the study in this direction is only on its 
very beginning. Many efforts should be made to set up a self-contained inflation theory 
on noncommutative spacetime. We still have little knowledge on thermodynamics of 
the deformed radiation field. 

From the linearity of the equations of electrodynamics on commutative spacetime, 
we know that photons do not interact with one another. The thermodynamics of 
the photonic field can be described well by an ideal Bose gas system. However, the 
photons in the deformed radiation field have strong correlations among themselves. 
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The strong correlations among the photons are caused by the nonlinearity of equations 
of electrodynamics on noncommutative spacetime. It is a difficult task to deal with 
the interactions among the photons. Numerous papers have been published [33]- [46] on 
the thermodynamics of the deformed photonic field. However, we still have not got a 
satisfactory distribution function of the deformed Bose gas. 

In this paper, we have treated the deformed photonic field as a strongly corre- 
lated system by making use of the Green's function theory. Equations of motion 
of the Green's function for the deformed photonic field was set up. The Tyablikov 
procedure [57] was used to decouple these equations of motion. A general form of sta- 
tistical distribution function for the deformed photonic field was got systematically. 
Based on this distribution function, we calculated the free energy, the entropy and the 
specific heat of the deformed photonic field in detail. The equation of state was ob- 
tained. We found that the Planck's formula for the deformed Bose gas is really different 
from the ordinary one. We got terms in Planck's formula, which correspond to the in- 
teractions among photons. Similar to the case of an ideal Bose gas, the total energy of 
the deformed Bose gas is proportional to the fourth power of the temperature. But the 
Stefan- Boltzmann's constant is effectively reduced by the noncommutative deformation 
of geometry. We noted that the density of the spectral frequency distribution of the 
energy of the deformed Bose gas has a maximum at a frequency cu m (/3ftu; m ~2.8), which 
is roughly independent of deformation parameter. Thus the displacement law is still 
remained for the deformed Bose gas. The results obtained in the paper can be used to 
investigate properties of the early universe (e.g. inflation model on noncommutative 
spacetime). All of these works are in progress. 
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Abstract 

By making use of the double-time Green function technique, we study ther- 
modynamics of a deformed Bose gas, which describes well properties of density 
intensive photonic gas and radiation fields of the early universe. General form of 
statistical distribution function is obtained. We show explicitly the expression 
of the distribution function in some limitation cases. The free energy, equation 
of state, specific heat and other thermodynamic properties of the deformed Bose 
gas are presented. 
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1 Introduction 

There are two main motivations to study a deformed many-body system. The first one 
comes from fundamental physics. It is well-known that the noncommutative geometry 
plays a very important role in unravelling the properties of the Planck scale physics. 
It has for a long time been suspected that the noncommutative spacetime might be 

a realistic picture of how spacetime behaves near the Planck scale[l]. In fact, the 
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noncommutative geometry naturally enters the theory of open string in a background 
-B-field[2]. In particular, the noncommutative geometry makes the holography[3] {e.g. 
the AdS/CFT correspondence) of a higher dimensional quantum system of gravity and 
lower dimensional theory possible. A kind of special regularization with exponentially 
increasing spacetime cutoff has been introduced based on the noncommutative geome- 
try. A very small displacement of the noncommutative deformation parameter from its 
classical value reduces sharply the entropy of quantum system of gravity. An adequately 
adopted noncommutative deformation of geometry makes the holography of a higher 
dimensional quantum system of gravity and lower dimensional theory possible. It was 
also discovered that simple limits of M theory and superstring theory lead directly to 
the noncommutative gauge field theory [4, 5]. The noncommutative field theory has 
been intensively studied in the past two decads[6, 7]. The most realistic laboratory for 
testing the Planck scale physics is the early universe and the black hole physics. From 
the standard picture of cosmology, the universe started from a very hot, dense phase 
about 15 billion years ago in the big bang. The very early universe was opaque due to 
the constant interchange of energy between matter and radiation. Several problems of 
standard big bang cosmology, such as the horizon, homogeneity, entropy and flatness 
problems, can be explained under the assumption that the very early universe under- 
went a period of extremely rapid inflation[8], driven by the potential of some assumed 
inflation fields. In particular, the inflation scenario is also able to explain the observed 
perturbations in the cosmic microwave background radiation [9, 10], namely as origi- 
nating ultimately from quantum fluctuations of the inflation field. However, in most of 
the current models of inflation, the physical length of perturbations was much smaller 
than the Planck length at the beginning of inflationfll, 12]. It is then of interest to 
investigate whether the predictions of the spectrum of cosmological perturbations are 
sensitive to the unknown super- Planck- scale physics[13]-[16]. Recent studies[17]-[24] 
by physicists and astrophysicists show that without strange inflation fields but just 
the noncommutative geometry and the deformed radiation dominated universe subject 
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to the noncommutative spacetime quantization indeed give rise to inflation as the ra- 
diation temperature exceeds the Planck temperature. Based on the noncommutative 
geometry and the deformed radiation, one can even get an alternative explanation for 
the horizon, homogeneity, entropy and flatness problems in the standard cosmology. 

Another motivation comes from quantum optics as well as molecular and atomic 
spectroscopy. A fascinating development of quantum optics in last two decades has 
taken place within two specific areas, one concerning the study of nonclassical states 
of light, and the other studying the interaction of light with matter. The nonclassical 
states of light offer the possibility of the so-called reduction of quantum noise below the 
standard limit, determined by the Heisenberg's uncertainty relation. Study on interac- 
tion of light with matter shows that the old Bohr quantum jumps are real. Dynamics of 
deformed photonic field interacting with matter presented many interesting properties 
both in the two important areas [2 5] -[29]. Great success has also been made by making 
use of a deformed oscillator to study molecular spectrum[30]-[32]. 

In this paper, we study thermodynamics of a deformed Bose gas, and in particular, 
of a deformed radiation field. On the topic, many papers have been published in more 
than ten years [33]- [46]. However, we still have not a general form of statistical distri- 
bution function of the deformed Bose gas, which should reduce to the Bose-Einstein 
distribution function for an ideal Bose gas in the q — 1 limit. In this situation, it is dif- 
ficult to investigate systematically the influence of the deformation on thermodynamic 
characteristics of a Bose system, which is interesting in fundamental physics and other 
subjects as discussed above. For this purpose, by making use of the double-time Green 
function technique, we obtain an algebra equation, which should be satisfied by the 
distribution function of a deformed Bose gas. The equation for distribution function 
can not be solved exactly in general. The situation is similar to the case of magnetic 
lattices. From this point of view, a deformed Bose gas is really a strongly correlated 
system. In the case of q — > 1, we give a useful expression of the distribution function 
for a deformed Bose gas. Thermodynamic properties, such as internal energy, specific 
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heat, free energy and equation of state, are discussed in detail. 

The paper is organized as follows. In section 2, a deformed Bose gas is discussed. 
We get the distribution function for a deformed Bose gas by using the double-time 
Green function technique in section 3. In section 4, we calculate some thermodynamic 
quantities for a deformed photon gas and compare them with those for the ideal Bose 
gas. Concluding remarks are given in section 5. 

2 A deformed Bose gas system 

The Lagrange density of a scalar field 4>(x) with mass m reads 

h 2 d(fr d<p 1 2 2 2 

£{x) = — — — m c <p . 1 

v ; 2 dx^ dx» 2 r K J 

The Euler-Lagrange equation immediately leads to the Klein-Gordon equation 



The canonically conjugate field is 



d(p(x) 



which leads to the Hamilton density 

H{x) = \ (^tt 2 (x) + (W0(x)) 2 + m 2 c 2 2 (x)) . (4) 
The Fourier decomposition of the scalar field is of the form, 

(j)(x)= [ l d3k ° (a(k)e- ikx + aUk)e ikx ) , 

= I $k ih 2 u k ( a{k)e - ikx + a \ {k y kx \ 
J - /(27r) 3 2^ fe V ' 



where we have used the notation cu k = £yj (hk) 2 + m 2 c 2 . The canonical quantization 
of the scalar field reads 

[0(x,t),7r(y,t)]=^c5 3 (x-y) . (6) 



In terms of the Fourier modes a(k) and at(k), the commutation relations can be rewrit- 
ten into the form 

[a(k),at(k')] = 5 3 (k - k') . (7) 
It is not difficult to express the Hamiltonian in terms of these Fourier modes, 

H = -J Sk%u k (a{k)a\k) + a\k)a(k)) . (8) 

In the quantum field theory, the vacuum state is defined as follows 

a(Jfe)|0) = 0. (9) 

And the iV-particle Fock space can be written down directly as 

WkJMh), n(k m )) = v ^ v > y ■ ■ ■ v ) = |o) • (io) 

\jn[k x )\ \jn[k 2 )\ \jn[k m )\ 

To measure how many particles of a certain momentum exist in the field system, one 
can introduce the particle number operator N(k), 

N{k) = a\k)a(k) . (11) 

It is easy to check that 

N(k i )\n{k 1 ),n(k 2 ), ■ ■ ■ ,n(k m )) = n(k i )\n(k 1 ),n(k 2 ), ■ ■ ■ ,n(k m )) . (12) 

Thus, a field can be seen as a many-body system that consists of free harmonic oscilla- 
tors with different momentum. Properties of the field system can be described exactly 
by the Weyl-Heisenberg algebra, 

[a, a 1 "] = 1 , 

(13) 

[N, a] = -a , [N, a f ] = a) , 

with the Fock space \n) = ^=r|0) as its representation. All properties , in particular, 
the thermodynamic properties of the many-body system can be investigated system- 
atically, by making use of the Hamiltonian (8) and the Weyl-Heisenberg algebra (13). 
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The natural quantum of the quantum field theory on noncommutative geometry 
is the so-called deformed harmonic oscillator. On the topic of deformed oscillators, a 
lot of papers have been published[47]-[49]. Similar to the case of an ordinary quantum 
field system, the many-body properties of a noncommutative quantum field system can 
be described by a collection of the deformed harmonic oscillators, 

H q = \ E huJ k (a q (k)al(k) + a\(k)a q (k)) . (14) 

The deformed harmonic oscillator is related with the simple harmonic oscillator as 
follows [50], 

S a , a; = atS . (15) 



a q = 



q~ — q 

where we have used the notation [x] q = r- 

q — q 1 

The deformed Weyl-Heisenberg algebra reads 

a q a\ - qa\a q = q~ N , 

W, aj] = 4 ' t^' ^q] = ~a q ■ 
The representation of the deformed Weyl-Heisenberg algebra is obtained by construct- 
ing the Fock space, 

yW (17) 
a q \n) q = ^J\n\ q \n - l) q , a\\n) q = yj[n + l] q \n + 1), , 

where [n] q \ = [n] q [n — l] q ■ ■ • [2] q [l] q . It is easy to check that one still can measure how 

many particles of a certain momentum exist in the deformed field system by making 

use of the particle number operator N(k), 

N(k i )\n(k 1 ),n(k 2 ), ■ ■ ■ ,n(k m )) q = n{k i )\n(k 1 ),n(k 2 ), ■ ■ ■ ,n(k m )) q . (18) 

The Hamiltonian of a deformed harmonic oscillator [51] -[54] can be written as 

H q = ^([N] q + [N + l] q ). (19) 

Below we will calculate the distribution function of a deformed harmonic oscillator 
system by making use of the double-time Green function technique. 
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3 Distribution function 

The double-time Green function (retarded or advanced) for a Bose system [55, 56] is 
generally defined as 

G p (t - o =~0p+ l)0(t - + (p - l)0(f - 0) ( W), fl(01> - (20) 

where we have used the notations, 

i(*) = exp(if)iexp(=ip) , 

(6) = Z- lr Tr(e-P H d) , Z = Tr(e^ H ) , 

here p = 1 for the retarded one and p = — 1 for the advanced one. 
In Eq.(20), setting 

A(t) = a(t), s(0 = at(0, p = l, 

(the retarded Green function is assumed), we get 

G R (t - t') = -i0(t - t')(a(t)a\t') - a\t')a{t)) . (21) 
The equation of motion of the retarded Green function reads 
ih^-G R (t-t') = h (8{t - t'){a(t)a\t') - a\t')a(t)) + 9(t - t') (^^^(f) - a\t') 



da(t) 



dt 
(22) 



The Heisenberg equation for the operator a(t) 



gives, 

= -if + <T^ +1) ) a(t) . (23) 

Then, we have 



ih^G R (t - H) = fiS(t - t') + (q^ + q~^) G R (t - t') . (24) 



Making use of the Fourier transformation 

G R (t - = ^G R {E)e-^-^ h , (25) 

we get 

ft 

= ~ Mg^+g-^ • (26) 
The spectrum theorem of Green functions tells us 

(a (t )a(t)) - i J ^ dE _ i} . (27) 

Then, we can write formally the distribution function of a deformed Bose gas as 

n(e) = IN) = — 77T^ • (28) 

exp( ^^> +g -^>) ) _ 1 

Evidently, in the case of q — >1, n(e) recovers the familiar Bose-Einstein distribution 
function. In principle, all interesting properties of a deformed Bose gas can be ob- 
tained by solving this equation for distribution function. This is to be one of the most 
advantageous features of the Green function approach. 

Now we try to solve Eq.(28) to gain an insight into the influence of the deformation 
on the distribution function and to compare the deformed distribution function with 
the Bose-Einstein distribution function. 

Expanding the right hand side of Eq.(28) in the neighborhood of q = 1 to the second 
order, we get an exactly solvable equation for the deformed distribution function. The 
solution is of the form, 

+ 



5 2 f3eet 



^- £ \J(5 2 fiee^ + ( e /fe - l) 2 ) 2 - 5 2 (3ee^ {5 2 {3ee^ - 2{eP e - 1)) , 

(29) 

where we have used the notation 5 = q — 1 for short. 

We note that it is difficult to compare the distribution function of the deformed 
Bose gas, Eq.(29), with the Bose-Einstein distribution function directly. For the aim, 
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we rewrite Eq.(29) as follows, 

(e* _ 1) _ Wee* (d 2 (3ee^ - 2 (e* - l))' 



n e 



(30) 



(e* - l)^ + F/fee* (a^/fee* + 2 (e/3e _ 1)2 y 

One can further expand the distribution function into a series over deformation pa- 



rameter 8, 



n(e) 



1 



-8 



, , I3ee li£ I3ee' i£ pee 



/3e 



(31) 



V 2 ^" 1 ) 2 (e^-1) 3 2(e^-l)V ' 
From the above form of the distribution function of the deformed Bose gas, we see 
clearly that, at the case of q = 1, the Bose-Einstein distribution is recovered. To 
further investigate the effect of deformation on the distribution function, we present a 
plot of the distribution function in Fig.l. 




500 
400 
300 
200 
100 

0.01 0.02 0.03 0.04" 

FIG.l The distribution function of a deformed Bose gas 

( start from up q = 1.000, 1.001, 1.002, 1.003, 1.004, 1.005, 1.006, 1.010, here we have used the notation 0s = E.) 

To describe behaviors of the distribution function of the deformed Bose gas near the 
ground state more accurately, we give another plot of the deformed distribution function 
in Fig.2. 



9 



8000 



6000 



4000 



2000 




o.ooi 0.002 0.003 0.004 FIG 2 T he behaviors of 
distribution function 
of a deformed Bose gas near the ground state 

(start from up q = 1.000, 1.001, 1.002, 1.003, 1.004, 1.005, 1.006, 1.010, here we have used the notation fie = E.) 



4 Thermodynamic properties 

In this section, we concentrate on the study of thermodynamic properties of a deformed 
photonic field. As being well-known, the dispersion relation of a photon is 



e = cp . 



(32) 



We know that a photonic field can be described by a collection of simple harmonic 
oscillators and have relation 



h = hu In + 



l 



One can then read out that 



hu ( 1 
P=-(n+- 



(33) 



(34) 



This relation and the Hamiltonian (19) give the dispersion relation of a deformed Bose 
gas, 



e = -hu 
2 



c 1 
h^ P+ 2 



+ 



hu P 2 



(35) 
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Expanding the right hand side of (35) in the neighborhood of q — 1 to the second 
order, we get 

S 2 5 2 , , 

£ = CP -24 CP+ W^ CP ■ (36) 
The density of energy levels is as follows, 

8nV SVWS + + 27Be 2 (-2 x 6^ A + (iSVBe + 2^3(4^ + 27 Be 2 )f") 



a(e) 



3h^ AB^AA* + 27W 2 ^ £ + 2y / 3{AA s + 27Be 2^ 1/3 

1/3 



9VBe + ^/3{AA 3 + 27Be 2 ) + &^A [18VBe + 2^3(AA 3 + 27Be 2 

X 2 

(9^^ + ^3(4^3 + 27Be 2 )) 

(37) 

where we have used the notations A = (1 -)c and .£> = — o-^-c 3 . 

V 24 ; 6ft V 

The mean value of the total energy of the system has the form, 



E 



/ en(e)a(e)de . (38) 
Jo 



To calculate E, we expand the right hand side of (37) to the third order in the neigh- 
borhood of q — 1 and rewrite the density of energy levels as the form 

Then we get 

where ((n) is the Riemann zeta function. 

We can further calculate the grand partition function In S through the thermody- 
namic relation E — — InS , 

op ' 

8n 5 V nVS 2 f8n 2 llvr 4 >/0 A /Im v, 

kS = " W (- + + «>C(3) J <^) 3 • (41) 

It is a straightforward calculation to obtain thermodynamic functions by making use 
of the grand partition function. The pressure of the deformed Bose gas is 
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The free energy of the deformed Bose gas is of the form, 



+ irV5 2 



six 



45 (he) 3 ' " ' " V 9 
The entropy of the deformed Bose gas is as 
U — F 



Uir 4 
~A5 



+ 



16C(3)1 



J (he) 3 



S 



3 ^ k4T3 v-^( 3 -^ + ^ + m s)) k ^v. 



(43) 



(44) 



T 45 h 3 c 3 V 9 45 7 

The first term in the above equation is just the entropy of an ordinary Bose gas. And 
the terms following it denote the effect of noncommutative deformation of geometry 
on the entropy. It is clear that the effect can not be ignored at least in the high 
temperature region, where the early universe lived. We give a plot of the entropy vs 
temperature for different deformation parameters in Fig. 3. 



5-10 



4 • 10 



3-10 



2 • 10 



1 • 10 




FIG. 3 The entropy per unit volume of a deformed Bose gas 

( start from up q = 1.0, 1.1, 1.2, 1.3, 1.4). 



Finally, the specific heat of the deformed Bose gas is 



32tt 5 V k A T 3 



'32n 2 Uir 4 
+ 



+ 192C(3) 



k 4 T 3 



(45) 



15 h 3 c 3 " ' " \ 3 ' 15 ' "~^"'J h 3 c 3 ' 
At the case of q — 1, we recover the result for an ordinary Bose gas. To show the effect 
of deformation on the specific heat, we give a plot of the specific heat per unit volume 
of a deformed Bose gas vs temperature for different deformation parameters in Fig. 4. 

From the above discusses, we can also obtain the Planck's formula for the deformed 
Bose gas, 
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FIG. 4 The specific heat per unit volume of a deformed Bose gas 

( start from up q = 1.0, 1.1, 1.2, 1.3, 1.4). 

hcu 3 l7S 2 hu 3 5 2 /mVe^ 



u(u,T) = 



u 



c 3 7r 2 (e^ - 1) 24c 3 7r 2 (e^ - 1) 2c 3 7r 2 (e^ - l) 2 

5 2 f3h 2 u A e^ 5 2 (3h 2 u 4 e^ 
c 3 7r 2 (e^ - l) 3 2c 3 7r 2 (e^ - l) 4 ' 




(46) 



1 . 4 
1.2 

1 

. 8 
. 6 
. 4 
0.2 

~ 4 6 8 10 " FIG. 5 The Planck's 

formula of a deformed Bose gas 

( start from up q = 1.0, 1.1, 1.2, 1.3, 1.4, here we have used the notation / = fihu>.) 

In Fig. 5, we give a plot of the Planck's formula for the deformed Bose gas. The first 
term in equation(46) describes an ideal gas consisting of photons. The linearity of 
the equations of electrodynamics on commutative spacetime guarantees that photons 
do not interact with one another. However, the other terms in this equation denote 
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interactions among photons or nonlinearity of equations of electrodynamics caused by 
noncommutative deformation of spacetime. At high frequencies {f3Twj ^> 1), we have 

hco 3 e-^ 175 2 hu; 3 e-P h " ^Ve"^ 

(47) 



u(uj, T) 



C 3 7T 2 



24c% 2 



2c% 2 



C 3 7T 2 2c 3 7T 2 

This is the Wien's formula of the deformed Bose gas. We give a plot of the Wien's 
formula for different deformation parameters in Fig. 6. In spite of a complex form of 
the Planck's formula and the Wien's formula, we note that the density of the spectral 
frequency distribution of the energy of the deformed Bose gas has a maximum at a 
frequency u m (f3hu m ~2.8), which is roughly independent of deformation parameter. 
Thus the displacement law is still remained for the deformed Bose gas. 




4 6 8 10 

FIG. 6 The Wien's formula of a deformed Bose gas 

( start from up q = 1.0, 1.1, 1.2, 1.3, 1.4, here we have used the notation / = (3hui). 



In the opposite limiting case of low frequencies ((3Tilo <C 1), we get the Rayleigh- Jeans 
formula for the deformed Bose gas, 



u(uj,T) = 



kTu 2 175 2 kTu 2 5 2 kTuj 2 ,hu 



2 24c% 2 

2(1 m\9 



2cH 2 y kT 



(77F + 1 ) 



5 2 {kT) 2 u huj_ ^ 5 2 {kTf (fiuj N 



(48) 



C 3 7T 2 /l ^kT 



2cH 2 h 2 \ kT 
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Fig. 7 shows a graph corresponding to the Rayleigh- Jeans formula with different defor- 
mation parameters. 

u 




f 



0.2 0.4 0.6 0.8 1 

FIG. 7 The Rayleigh-, Jeans formula of a deformed Bose gas 

( start from up q = 1.00, 1.06, 1.12, 1.18). 

We obtain the Stefan- Boltzmann's law for the deformed Bose gas of the form, 




0.2 0.4 0.6 0.8 1 



FIG. 8 The Stefan- Boltzmann law of a deformed Bose gas 

( start from up q = 1.0, 1.1, 1.2, 1.3, 1.4). 

Similar to the case of an ideal Bose gas, the total energy of the deformed Bose gas 
is proportional to the fourth power of the temperature. But the Stefan- Boltzmann's 
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constant is effectively reduced by the noncommutative deformation of geometry. The 
fact gives also a signal that the deformed Bose gas is a strongly correlated system. We 
give a plot of the Stefan- Boltzmann's law for different deformation parameters in Fig.8. 

In an adiabatic expansion (or compression) process, the entropy of the system is a 
constant and the valume and temperature are related by VT 3 =constant. From (42), 
we can get the equation of state, 

p\/ 4/3 = c q . (50) 

5 Concluding Remarks 

The noncommutative geometry plays a very important role in the study of the Planck 
scale physics. Quantum field theory and superstring theory (or M theory) on noncom- 
mutative spacetime have been studied for a long time. It is well-known that the most 
realistic laboratory for testing the Planck scale physics is the early universe. The very 
early universe was opaque due to constant interchange of energy between matter and 
radiation. It is clear that the radiation field at the early universe can not be described 
well by an ideal Bose gas system. There should be strong correlations among photons 
at this stage of evolution of the universe. In fact, attention has been paid recently to 
the investigation on the deformed radiation dominated universe. Preliminary results 
show that spacetime quantization gives rise to inflation as the radiation temperature 
exceeds the Planck temperature. Of course, the study in this direction is only on its 
very beginning. Many efforts should be made to set up a self-contained inflation theory 
on noncommutative spacetime. We still have little knowledge on thermodynamics of 
the deformed radiation field. 

From the linearity of the equations of electrodynamics on commutative spacetime, 
we know that photons do not interact with one another. The thermodynamics of 
the photonic field can be described well by an ideal Bose gas system. However, the 
photons in the deformed radiation fields have strong correlations among themselves. 
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The strong correlations among the photons are caused by the nonlinearity of equations 
of electrodynamics on noncommutative spacetime. It is a difficult task to deal with 
the interactions among the photons. Numerous papers have been published [33]- [46] on 
the thermodynamics of the deformed photonic field. However, we still have not got a 
satisfactory distribution function of the deformed Bose gas. 

In this paper, we have treated the deformed photonic field as a strongly corre- 
lated system by making use of the Green's function theory. Equations of motion 
of the Green's function for the deformed photonic field was set up. The Tyablikov 
procedure [57] was used to decouple these equations of motion. A general form of sta- 
tistical distribution function for the deformed photonic field was got systematically. 
Based on this distribution function, we calculated the free energy, the entropy and the 
specific heat of the deformed photonic field in detail. The equation of state was ob- 
tained. We found that the Planck's formula for the deformed Bose gas is really different 
from the ordinary one. We got terms in Planck's formula, which correspond to the in- 
teractions among photons. Similar to the case of an ideal Bose gas, the total energy of 
the deformed Bose gas is proportional to the fourth power of the temperature. But the 
Stefan- Boltzmann's constant is effectively reduced by the noncommutative deformation 
of geometry. We noted that the density of the spectral frequency distribution of the 
energy of the deformed Bose gas has a maximum at a frequency cu m (/3ftu; m ~2.8), which 
is roughly independent of deformation parameter. Thus the displacement law is still 
remained for the deformed Bose gas. The results obtained in the paper can be used to 
investigate properties of the early universe (e.g. inflation model on noncommutative 
spacetime). All of these works are in progress. 
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